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The upper mass bound of the lightest neutral Higgs scalar is studied in the \i problem 
solvable extra U(l) models by using the analysis of the renormalization group equations. 
In order to restrict the parameter space we take account of a condition of the radiative 
symmetry breaking and some phenomenological constraints. We compare the bound ob- 
tained based on this restricted parameter space with the one of the next to the minimal 
supersymmetric standard model (NMSSM). Features of the scalar potential and renor- 
malization group equations of the Yukawa couplings among Higgs chiral supermultiplets 
are rather different between them. They can reflect in this bound. 
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1. Introduction 

Low energy supersymmetry is one of the main subjects of present particle physics. It 
is considered to solve a weak scale stability problem called the gauge hierarchy problem in 
the standard model (SM). Although we donot have any direct evidence for it, it has been 
stressed that the gauge coupling unification occuring in a rather precise way in the minimal 
supersymmetric extension (MSSM) of the SM may be an encouraging sign for the presence 
of the low energy supersymmetry. In the MSSM its phenomenology crucially depends on 
the soft supersymmetry breaking parameters and then it seems to be difficult to make 
useful predictions unless we know how the supersymmetry breaks down. However, there 
is an important exception that the lightest neutral Higgs scalar mass cannot be so heavy 
and it is mainly controled by the feature of the weak scale symmetry breaking fl[|. This 
is not heavily dependent on the feature of the soft supersymmetry breaking parameters 
at least at the tree level. Thus the knowledge of its possible upper bound is crucial to 
judge the validity of the low energy supersymmetry from a viewpoint of the energy front 
of the accelerator experiment. This aspect has been extensively studied taking account 
of a radiative correction mainly due to a large top Yukawa coupling 0. 

It is well known that there still remains a hierarchy problem called // problem in the 
MSSM. Why a supersymmetric Higgsino mixing term parametrized by fi is a weak scale 
cannot be explained in the MSSM 0. A simple and promising candidate for its solution 
is an extension of the MSSM by the introduction of an extra U(l) gauge symmetry and 
a SM singlet field S with a nonzero charge of this extra U(l) @, [5[]. The essential feature 
of this model is described by the following superpotential 

W u[iy = \SH X H 2 + kSgg + h t QH 2 T + • • ■ , (1) 

where H x and H 2 are usual doublet Higgs chiral superfields and the ellipses stand for the 
remaining terms in the MSSM superpotential other than the \i term and the top Yukawa 
coupling. In the second term g and g stand for the extra color triplet chiral superfields 
which are important to induce the /i scale. In the superpotential Wu(\y we also explicitly 
write the top Yukawa coupling because of its importance in the electroweak radiative 
symmetry breaking as in the case of the MSSM ||. 

The vacuum of these models is parametrized by the vacuum expectation values (VEVs) 
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H 2 ) = , (5) = u, (2) 



of Higgs scalar fields such as 

/ 

(Hi) = 



where v\ and v 2 are assumed to be positive and vf + v\ — v 2 {= (174GeV) 2 ) should be 
satisfied^]. The vacuum in this model is parametrized by tan/3 = t^/fi and u. The extra 
U(l) symmetry is assumed to be broken at the region not far from the weak scale by a 
VEV of the scalar component of S because of the radiative effect caused by the second 
term in Wuny @, |5| and then the \x scale is induced as \x = Xu. Thus in this model the 
sign of n is fixed as the one of u automatically. 

This extra U(l) symmetry forbids a bare \x term in the superpotential and simultane- 
ously makes the model free from the massless axion and tadpole problems. These features 
seem to make this model more promising than the next to the minimal supersymmetric 
standard model (NMSSM) 0] which is similar to this extra U(l) model but is extended 
only by a SM singlet chiral superfield S with the superpotential 

Wnmssm = \SH X H 2 + ^kS 3 + h t QH 2 T + ■■■. (3) 

It is also interesting that this kind of extra U(l) models can be often obtained as the 
effective models of a lot of superstring models |§. Various interesting features of this type 
of models have been studied in many works by now [9—12]. Among the phenomenology 
of these models the lightest neutral Higgs scalar mass is also an important target for the 
detailed investigation. Of course, also in these models the lightest neutral Higgs scalar can 
be expected to be generally not so heavy. The interesting point is that its upper bound 
can be calculable with no dependence on the soft supersymmetry breaking parameters at 
least at the tree level as in the case of the NMSSM [13—17]. The dependence on the soft 
supersymmetry breaking parameters comes in through the loop correction mainly due to 
the large top Yukawa coupling and the second term of Eq. (1). 

In this paper we estimate the upper bound m^o of this lightest neutral Higgs scalar 
mass on the correct vacuum. The correct vacuum is determined as the radiatively induced 
minimum of the effective potential in the suitable parameter space. In this approach we 
1 In the following discussion we donot consider the spontaneous CP violation. Under this assumption 
the sign of u cannot be fixed freely but it should be dynamically determined by finding the potential 
minimum. 
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use the one-loop effective potential and solve the relevant renormalization group equations 
(RGEs) numerically. We will pay our attention on the comparison of this upper bound 
with the one of the NMSSM within the phenomenologically allowable parameter region. 
In the NMSSM it has been known through many works that the triviality bound of a 
Yukawa coupling A of Higgs chiral superfields strictly control the upper bound of the 
lightest neutral Higgs mass |16| . Our approach is somehow different from this usual 
one. We find the phenomenologically acceptable parameter subspace in the rather wide 
parameter space by taking account of the radiative symmetry breaking condition and 
some phenomenological conditions such as the chargino mass and the charged Higgs scalar 
mass etc.. The estimation of the upper mass bound of the lightest neutral Higgs scalar is 
carried out in this restricted parameter subspace. Although the result of this approach is 
necessarily dependent on the assumption for the soft supersymmetry breaking parameters, 
we consider that it is possible to obtain the useful results by studying the wide region of 
the parameter space. 

2. Extra U(l) models 

In this section we discuss more detailed features of the extra U(l) models and give the 
basis of the present study. Since the NMSSM is well known and discussed in many papersf], 
it is cinvenient to explain the points by using the extra U(l) models. The superpotential 
of our considering extra U(l) models is defined by Eq. (1). Soft supersymmetry breaking 
parameters are introduced as 

+ (A x \SHiH 2 + A k kSgg + A t h t QH 2 T + h.c.) , (4) 



where the first two terms are mass terms of the scalar component <pi of each chiral super- 
multiplet and of gauginos A a . We use the same notation for the scalar component as the 
one of the chiral superfield to represent the trilinear scalar couplings in the last parenthe- 
ses. Other freedoms remaining in the models are extra matter contents and a type of extra 
U(l). On these points we confine our study into the typical extra U(l) models derived 
from Eg, which are listed in Table 1. At the TeV region they are assumed to have only one 
extra U(l) symmetry which is broken only by the VEV of S and give a solution to the 
problem ||11|| . As discussed in Ref. [16] for the case of NMSSM, the extra matter contents 
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Table 1 The charge assignment of extra U(l)s which are derived from Eq p2[ . These charges are 
normalized as = 20. 
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affect indirectly the low energy value of the Yukawa coupling A through the influence on 
the running of the top Yukawa coupling. This is rather important to estimate the Higgs 
mass bound. In the present model such kind of effects on the Yukawa couplings may also 
be expected but its effect is more complicated than the NMSSM as discussed later. If 
we introduce the extra field contents arbitrarily, the cancellation of the gauge anomaly 
may require to introduce the additional fields which again affect the running of Yukawa 
coupling A and so on. Thus for the estimation of the Higgs mass bound it is important 
to fix the matter contents in the anomaly free way in the present study. 

As the matter contents we assume the MSSM contents and additional extra matter 
fields 



3{Q,U,D,L,E) + (H 1 ,H 2 



MSSM 



3(g,g)+2(H 1 ,H 2 )+3(S) + 3(N), 



which can be derived from three 27s of Eq shown in Table 1. This set satisfies the anomaly 
free conditions. We can also add extra fields to these in the form of vector representations 
constructed from the fields listed in Table 1. Here we consider the following two cases as 
the additional extra chiral superfields 



(A) (H a ) + (H*) 



(B) (g + H a + H b ) + (g* + H: + HZ), 



where a, b = 1 or 2 and the fields in the second parentheses come from 27* of E e . At least 
on the sector of SU(3)c , xSU(2) i xU(l)y these matter contents are the same as the one 
of [MSSM + n(5 + 5*)] where 5 and 5* are the representations of the usual SU(5). The 
case (A) corresponds to n = 3 and (B) to n = 4. The n = 3 is the critical value for the 
one-loop /3-function of SU(3). It makes this one-loop /3-function be zero. The interesting 
point of these field contents is that the unification scale of SU(3)(7xSU(2) i xU(l)y is not 
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shifted from the MSSM one. The n = 4 case saturates the /3-function for the pertubative 
running of gauge couplings up to the unification scale ~ 3 x 10 16 GeV. Although this 
addition seems to be artificial, this type of spectrum can be expected in the Wilson line 
breaking scenario of the E 6 type superstring model. We use these contents to compare 
the feature between the NMSSM and our extra U(l) models^. 

The existence of multi-generation extra fields brings an ambiguity in Eq. (1). The 
coupling A and k can have generation indices for extra fields such as S, Hi, H 2 , g and 
g. On this point we make the following assumption to make the argument simple. f] Only 
one S can have the couplings in Eq. (1) and one pair of (Hi,H 2 ) corresponding to the 
one of the MSSM alone gets the VEVs. Extra colored singlets (gi,gi) have a diagonal 
coupling to this S as kiSg^i, where all the coupling constant k{ show the same behavior 
in the RGEs because (gi,g~i) are completely symmetric for the generation index i in the 
models. The fermion components of gi and gi can get mass through this coupling.^ On 
the other hand, the fermion components of the remaining S which donot couple to the 
usual Higgsinos in Hi and H 2 can get their masses through the one-loop correction. From 
a viewpoint of the model construction, the serious phenomenological problem will be 
how the fermion components in other remaining extra matter fields can get their masses. 
Although they can be generally massive through the gaugino mediated one-loop diagrams, 
their magnitude seems not to be enough to satisfy the phenomenological constraints. In 
the £_ model given in Table 1 we can introduce the intermediate scale through the D- 
flat direction of iV and iV* whose exsistence does not affect our discussion in the later 
part of this paper. If this is the case, they can have the weak scale masses through the 
nonrenormalizable interactions in the superpotential such as -j^j—NN*gg*. Although this 
is phenomeno logically important, it can be improved by the suitable extension without 
changing the following results and thus we donot get involved in this point further here. 

In our considering models the tree level scalar potential including the soft supersym- 

2 We should note that if every extra U(l) is broken near the unification scale, these models are equal 

to the NMSSM with the equivalent extra matters as discussed in Ref. p| . 

3 As far as we use Eq. (7) for the upper bound of the lightest neutral Higgs mass, this assumption seems 

to be reasonable. This assumption affects the RGEs of some parameters and also one-loop correction to 

the Higgs scalar mass. Other cases will be discussed later. 

4 If we change the charge assignment for some fields, some extra fields can be heavy at the intermediate 

scale as discussed in p7[ . Although this kind of possibility can be realized in the £± model, we donot 
consider it in this paper. 
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metry breaking terms can be written as 

V = ^(g 2 2 +g 2 )(\Hi\ 2 -\H 2 \ 2 ) 2 +(\XSHi\* + \XSH 2 \ 2 
+ mllH^ 2 + m 2 2 \H 2 \ 2 - {A X XSH X H 2 + h.c.) 

+ ~g% (Q^H.l 2 + Q 2 \H 2 \ 2 + Q S \S\ 2 ) 2 + X 2 \H X H 2 \ 2 + m 2 s \S\ 2 , (5) 

where Qi, Q 2 and Q$ are the extra U(l) charges of Hi, H 2 and S, respectively. The first 
two lines are found to have the corresponding terms in the MSSM if we remind the fact 
that /i is realized as fi = Xu. The third line contains new ingredients. Its first term is a 
D-term contribution of the extra U(l) and qe stands for its gauge coupling constant. 
Potential minimum condition for Eq. (5) can be written as, 

ml = —Agl + 9l){vl - v\) - \g 2 E Qi{Qivl + Q 2 v\ + Q s u 2 ) - A 2 (u 2 + v 2 2 ) + XA x u-, 
4 4 vi 

m l = -Ml + 9l)( v l - v l) - - A g%Q2(Qivl + Q 2 v\ + Q s u 2 ) - X 2 {u 2 + vl) + XA\u—, 
4 4 v 2 

ml = -]glQs(Qiv 2 + Q 2 v 2 2 + Q s u 2 ) - X 2 {v\ + v 2 2 ) + XA X ^. (6) 
4 u 

This constrains the soft SUSY breaking masses of Higgs scalars around the weak scale. 
As the second derivative of Vo in Eq. (5) we can derive the mass matrix of the CP-even 
neutral Higgs scalar sector which is composed of three neutral components H®, H® and 
S. The goodness of this treatment has been discussed in the MSSM case [[| and we 
follow this argument. If we note the fact that the smallest eigenvalue of any matrix is 
always smaller than any diagonal elements, we can obtain the tree level upper bound 
of this lightest Higgs scalar in an independent way of the soft supersymmetry breaking 
parameters by transforming the basis into the suitable one. This upper bound can be 
written as PI 1 



(0)2 ^ 2 

m h0 < m 5 



2A 2 a 2 
cos 2 23 + — ^ sin 2 2(5 + yE . f Qi cos 2 3 + Q 2 sin 2 

9i + 92 9i +gi K 



(7) 



where we used the potential minimization condition (6). The first two terms correspond 
to the ones of the NMSSM in which their behavior has been studied in many works 
rj, [T5|, [T(J. The running of a coupling constant A and its triviality bound have been 



shown to be crucially dependent on the extra matters |16|]. The extra U(l) effect appears 
through the last term which is its D-term contribution. Equation (7) can show the different 
tan/3 dependence from the one in the MSSM depending on the value of A and also the 



type of extra U(l). In the case of MSSM the upper bound of the lightest neutral Higgs 
mass always increases with tan/5 in the region of tan/? > 1. If A ^ 0.6, the present 
model also shows the same behavior. On the other hand, for the same region of tan /3 its 
upper bound can decrease with increasing tan/? when A > 0.6 which does not depend on 
the model so heavily. The NMSSM shows the similar feature, which can be seen in Ref. 
[I6| . Although this may be potentially altered by the radiative correction, it is one of the 
typical features coming from the XSH1H2 in these models different from the MSSM. 

We should note that the bound formula Eq. (7) is applicable only in the case of 
u ^> v 1, t>2 0. In the extra U(l) model the value of u can be constrained from below by the 
conditions on the mass of this extra U(l) gauge boson and its mixing with ordinary Z°. 
As far as we donot consider the special situation such as tan 2 (3 ~ Q1/Q2 under which the 
mixing with the ordinary Z° is negligible, the hierarchical condition u> Vi,V2 should be 
imposed to satisfy the phenomenological constraints on the extra Z' mass and its mixing 
with ordinary Z° [[L2|]. In the sufficiently large u case A may be constrained into a limited 
range required by the successful radiative symmetry breaking at the weak scale so that 
Xu(= n) takes a suitable value. In the NMSSM this kind of constraint on A is expected 
to be weaker than the one of the extra U(l) model since u has no phenomenological 
constraint at this stage. Anaway, we need the RGE study to check whether A can be 
constrained in a substantial way by this condition. 

3. The comparison of extra U(l) models and the NMSSM 

It is useful to discuss some qualitative features of the extra U(l) models and the 
NMSSM in more detail before comparing the mass bound of the lightest neutral Higgs in 
both models. Although the extra U(l) models and the NMSSM have the similar feature 
related to the /x term, they are expected to show rather different behavior in the running 
of Yukawa couplings k, k and A. The top Yukawa coupling has the same one-loop RGE 
in both models as, 

-P- = -\ Uhl + A 2 - ™g>) . (8) 
din// 16vr 2 V * 3 y V V ; 

In the present field contents g% takes larger value at Mx than the one of the MSSM. Even 
if the initial value of h t takes the large value like 0(1), the /3-function in Eq. (8) can be 
5 In the case of u < v\ , f 2 the diagonal element corresponding to S can be smaller than the right-hand 
side of Eq. (7). In such a case we cannot use Eq. (7) as the bound of the lightest neutral Higgs mass. 
We will exclude it from our study. 



8 



small due to the cancellation between a h t term and a g% term. As a result, h t tends to 
stay near its initial value at the intermediate scale independently whether it starts from 
a large value or a small value. This feature is shared by both models. On the other hand, 
the one-loop RGEs of k, k and A are largely different from each other. They are witten 
as, in the NMSSM, 

dn 



^( &2+6A '' 



<iln/i 167T 5 

dX X (3h 2 t +2K 2 + 4X 2 ). (9) 



dlnfi 167T 

and in the extra U(l) model, 

dk k „ s , n 16 



din, 16,U (3iV » + ^ + ^ ~ 
'' A A 'ih 2 , + 3/V a /fc 2 + 4A 2 ) , (.10) 



cHn/i 167T 2 

where N g is a number of the pair of the singlet colored fields g and g which have a coupling 
to S. In these RGEs we neglect the effect of gauge couplings 5% 9i and qe Q. At first we 
consider the running behavior of k and k. Since fc has an effect of #3, it can be rather 
larger at the intermediate scale than k which has no such effect and rapidly decreases 
according to lowering energy. This is important to determine the value of u realized in 
both models, which are mainly determined by m| at the low energy region. They are 
controled by the one-loop RGE as 



dm 2 s 



\ (2K 2 (3m 2 s + A\) + 2A 2 (m| + m 2 Hl + m 2 H2 + A 2 X )) (11) 



din 11 87r 2 
in the NMSSM and 

^ = g^a (?N g k 2 {m 2 s +ml + m 2 s + A 2 ) + 2A 2 (m| + m 2 Hl + m 2 H2 + A 2 X )) , (12) 

in the extra U(l) models. The larger k compared with k makes m| much more negative 
in the extra U(l) models. The larger value of u is expected in the extra U(l) models if 
we remind Eq. (6). 

As easily seen from the RGE of A in the extra U(l) model, the running of A is made 
fast by the existence of the second term of Eq. (1) which is needed for the successful 
6 In these equations we cannot find the fixed ratio point other than k — or A = as far as N g ^ 
even if we ignore g$. This is very different situation from the NMSSM which has been discussed in Ref. 
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radiative symmetry breaking of these models |4| . The one- loop /3-function of the coupling 
k has a contribution of g% differently from the case of k in the NMSSM. If we start k 
and k from the large values at the unification scale, this feature can keep k rather large 
at the intermediate region and then the running of A can be made fast by its effect 
compared with the one of k in the NMSSM. This feature tends to make the value of A 
at the low energy scale smaller compared with the NMSSM case if the same initial value 
is adopted at least. However, the initial value of k and k should be controled from the 
requirement of the radiative symmetry breaking from our view point since they play an 
important role in this phenomenon. We need the numerical analysis to study this aspect 
in more quantitative way. The extra matter effects on the RGEs are also rather different 
between the NMSSM and the extra U(l) models. As far as all the couplings are within the 
perturbative regime, the larger number of extra matter fields make the gauge couplings at 



the unification scale larger. As pointed out in fL6[ , in the NMSSM this indirectly makes 
the low energy value of A larger through the smallness of h t at the intermediate scale 
whose /3-function in Eq. (8) is kept small there. On the other hand, in the extra U(l) 
models the runnings of k and A are simultaneously affected by the extra matters in both 
direct and indirect manner, as is easily seen in Eq. (10). 

We know from these considerations that the resulting low energy values of A and u 
are rather different in both models. We should note that these values affect the upper 
bound of the lightest neutral Higgs scalar mass. Although Eq. (7) shows A is crucial to 
determine the tree level bound, u is essential to determine the magnitude of the one-loop 
effect, especially in the extra U(l) models. The radiative correction to Eq. (7) can be 
taken into account based on the one-loop effective potential. It is well-known that the 



one-loop contribution to the effective potential can be written as |L9|, ^0 



where Ai 2 is a matrix of the squared mass of the fields contributing to the one-loop 
correction and A is a renormalization point. In the usual estimation of the lightest neutral 
Higgs mass in the NMSSM the top and stop contributions to V\ are mainly considered as 
the relevant fields because of their large Yukawa coupling. However, in the study of the 
extra U(l) models k is rather large and then we should also take account of the effect on 
M 2 from the extra singelt colored chiral superfields g and g which have a coupling with 
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S. A mass matrix of the stops is written as 

hii)n(— At 4- A?/, rot ff\ \ 

(14) 



1 m 2 Q + lr^v\ h t v 2 (-A t + Am cot (3) 



h t V2(—A t + Xu cot (3) rfif + hfv 

and the one of the s-gquarks is expressed as 



m 2 g + k 2 u 2 —Akku + \kv\V2 
\ —Aj,ku + \kv1V2 Trig + k 2 u 2 



(15) 



where fh 2 Q f, m 2 g and A t , Ak are soft supersymmetry breaking parameters. Here a D-term 
contribution is neglected as it has been done in many previous investigations of the MSSM 
0. Mass eigenvalues of these mass matrices are respectively expressed as, 



™| = ^1 + ™t) + h t v 2 ± J\( rh Q- ™f) + hM(~A + Am cot f3) 2 , 



m 3i = ^( m 9 + m i) + ^ 2 "" 2 ^ y ^( m 9 _ m i) + (~A k ku + \kviv 2 ) 2 . (16) 

If we estimate the upper bound of the lightest Higgs mass in the same procedure as the 
one used to obtain Eq. (7) by minimizing the one- loop effective potential V e s —V0 + V1, 
the following one- loop correction should be added to the right-hand side of Eq. (7): 

Amh ° = 2 [M " C ° S P+ 2d^2~ Sm2(3+ 2 " Sm ^ (17) 

From these we find that u can crucially affect to the mass bound through the one-loop 
effect of gquark sector in the extra U(l) models. This addtional effect cannot be escapable 
as far as the occurence of the radiative symmetry breaking is required. 

It may also be important to take account of the difference in both models coming from 
some phenomenological constraints, in particular, the ones related to A and u. Although 
this kind of constraints depend on the values of soft supersymmetry breaking parameters, 
it may be useful to improve the upper bound estimation based on the triviality bound of 
A. We should remind the fact that the chargino mass, the charged Higgs mass and squark 



masses are dependent on A and u ||11||. The chargino and the charged Higgs scalar have 



the same constituents as the MSSM. However, they have a different mass formulas from 



the MSSM ||12|| . In both models the chargino mass is expressed as 
m x ± = - (A V + 2m 2 v + M 2 



±J- (2m 2 w cos 2(3 + \ 2 u 2 - M$f + 2m 2 w (-Aw sin (3 + M 2 cos pf 
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where raw and M 2 represent the W boson and the gaugino A^ masses. The charged 
Higgs scalar mass has the different mass formula between both models. In the extra U(l) 
models it is expressed as 

2A 2 \ 2A x \u 
Of) + sin 2/? 



m 2 H± = m 2 w ( 1 - — ] + — X —, (19) 



while in the NMSSM it is written as 

, 2 / 2A 2 \ 2(A x \u-k\u 2 ) , . 

Recently the lower bounds of these masses become larger and we may use these to put 
some constraints on A and u. Another important point to use Eq. (7) is that it must 
be smaller than other two diagonal mass matrix elements of the 3x3 neutral Higgs 
scalar mass matrix. Especially the diagonal mass for the singlet Higgs scalar S can give 
a substantial constraint on u. Its tree level formula is 

2 1 2^,2 2 , A X \VlV 2 / 01 x 

m m = ~9eQs u + (21) 

Hi 

in the extra U(l) models, while it is expressed as 

9 2 2 A\\VlVo 

m 2 o = 4kV + — - — — - A k ku 22 

3 U 

in the NMSSM. This constraint may be substantial in the NMSSM where there is no other 
clear constraint on the small u. 

4. Numerical analysis and its results 

In this section we numerically estimate the bound of m 2 i0 (= rn^ 2 + Am| ) by solving 
the RGEs and taking account of the phenomenological constraints presented above. In 
order to improve the one-loop effective potential |2lJ We use two-loop RGEs for dimension- 
less coupling constants and one-loop ones for dimensional SUSY breaking parameters, for 
simplicity. In this estimation we adopt the following procedure. As the initial conditions 
for the SUSY breaking parameters we take 

m\ = ( 7i m) 2 , M a — M, A t = A k = A K = A x = A, (23) 

where m 2 is the universal soft scalar mass and we introduce the nonuniversality repre- 
sented by 7j only among soft scalar masses of Hi, H 2 and S. We comment on this point 
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later. These initial conditions are assumed to be applied at the scale where the coupling 
unification of SXJ(2) L and U(l)y occurs. We donot require the regolous coupling unifi- 
cation of SU(3)c but only impose the realization of the low energy experimental value 
following Ref. (T^|. For the extra U(l) coupling g% we use the same initial value as the one 
of U(l)y at the unification scale M X - The initial values of these parameters are surveyed 
through the following region, 

< ht < 1.2 (0.1), < k, \k\ < 2.0 (0.2), < A < 3.0 (0.2), 

< M/M s < 0.8 (0.2), < m/M s , \A\/M S < 3.0 (0.3), (24) 

where in the parentheses we give the interval which we use in the survey of these parameter 
regions. Since the sign of k and A affect the scalar potential, we need to investigate both 
sign of them. We also assume that the RGEs of the model are changed from the ones of the 
supersymmetric extra U(l) models to the nonsupersymmetric ones at a supersymmetry 
breaking scale Ms for which we take Ms = 1 TeV as a typical numerical value |14|, |n§B 



As a criterion for the choice of the correct vacuum, we impose that the radiative 
symmetry breaking occurs correctly. We check whether the potential minimum satisfying 
the conditions such as Eq. (6) improved by the one-loop effective potential can satisfy the 
phenomenologically required conditions such as v = 174 GeV and m t = 174 GeV starting 
from the above mentioned initial conditions. It is not so easy to find this solution under 
the completely universal soft breaking parameters so that in our RGEs analysis we allow 
the nonuniversality in the region 0.8 < 7« < 1.2 among soft supersymmetry breaking 
masses of Higgs scalars. The nonuniversality of soft scalar masses are generally expected 
in the superstring models [p2]j . This treatment seems to be good enough for our purpose 
such as to estimate the upper mass bound of Higgs scalar. We also additionally impose 
the following phenomenological conditions. 

(i) m^ should be smaller than other diagonal components of the Higgs mass matrix (see 
also footnote 3 and the discussion related to Eqs. (13) and (14)). 

(ii) the experimental mass bounds on the charged Higgs bosons, charginos, stops, gluinos 
7 In principle wc should solve the RGEs of soft supersymmetry breaking parameters under the initial 

values given in Eq. (15) in order to estimate this scale Ms- However, we donot take such a way here, for 
simplicity. It is beyond the present scope to study the dependence of our results on the supersymmetry 
breaking scale Ms- 
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Fig.l Scatter plots of the radiative symmetry breaking solutions in the (tan /?, u) plane for the NMSSM 
1(a) and the £_ model 1(b). Solutions for the different h t (Mx) are classified. 

and Z' should be satisfied. Here we require the following values: 



m H ± > 67 GeV, m x ± > 72 GeV, rh tl 2 > 67 GeV, 
M 3 > 173 GeV, m z , > 500 GeV. 



(25) 



(iii) the vacuum should be a color conserving one p3| . 

We adopt only the parameters set satisfying these criterions as the candidates of the 
correct vacua and calculate the Higgs mass bound m^ for them. 

At first in order to see the difference in the allowed vacuum between the NMSSM and 
the extra U(l) models we plot the radiative symmetry breaking solutions for the present 
parameter settings in the {tan ft, u) plane in Fig. 1. Solutions are classified by the initial 
value of ht at Mx into three classes which show rather different qualitative features. As 
an example of the extra U(l) models we take the £_ model here but the t] model has 
been checked to show the similar feature to the £_ model. We take the case (A) as the 
extra matter contents. Throught the present calculation an effect of the translation of 
the running mass to the pole mass p4| is taken into account to determine tan (3. We take 
tan/5 < 15 and neglect the large tan/5 solutions since the bottom Yukawa coupling is 
assumed to be small in the RGEs so that in the present analysis the large tan ft solutions 
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Fig. 2 Scatter plots of the radiative symmetry breaking solutions for the NMSSM and the £_ model in 
the (tan/3, k or k) plane 2(a) and the (tan/3, A) plane 2(b). The values of k, k and A are the ones at mt- 

cannot be recognized as the appropriate ones. Figure 1 shows that the £_ model can have 
solutions in the larger u region of the (tan /3, u) plane compared with the NMSSM. As 
mentioned in the previous section, this is a result that k can be larger than k at the m t 
scale due to the SU(3) effect. This is shown in Fig.2, where the values of k(m t ), K(m t ) 
and A(mt) corresponding to each solution are plotted for tan/3. The soft scalar mass m| 
of the singlet Higgs scalar S becomes much more negative in the extra U(l) models than 
in the NMSSM. In the sufficiently large u region the potential minimum condition for u 
reduces to 



4m< 



u 



9lQ 2 s 



for extra U(l), 



mt 



u 



2k 1 



for NMSSM. 



(26) 



In the NMSSM u depends not only on m| but also on k and as a result u can take a 
rather large value. In the £_ model the smaller u region such as u < 1 TeV is cut due 
to the experimental extra Z mass bound. Also in the NMSSM very small u seems to be 
forbidden. This seems to be a result of the phenomenological conditions (i) and (ii). 

The big qualitative difference of the vacuum in both models is that there can be large 
u solutions for tan/3 ^ 5 corresponding to ht(Mx) = 0.3 in the extra U(l) model. One 
reason of this is that the smaller A(m t ) is realized in the extra U(l) models than in the 
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Fig. 3 Boundary values of m h o of the lightest neutral Higgs mass as a function of tan (3 in the NMSSM. 
Full data are used to draw 3(a). In 3(b) we impose rh(Mx) = 1 TeV. All solutions satisfying 2.4 TeV 
< u < 2.6 TeV are also plotted in 3(b). 

NMSSM. This is clearly shown in Fig. 2(b). The discussion on this aspect has already 
given based on the RGE in the previous section. On this point we should also note 
that in the tan/5 > 5 region the small X(mt) is allowed. Thus \l — \u can be in the 
suitable range even if u is large. However, the boundary value of u seems not to have 
so strong dependence on \{mt) in both models and the value of \u does not seem to 
be strictly restricted by the radiative symmetry breaking at least within the parameter 
region searched in this paper. 

In Figs. 3~5 we give the results of our numerical estimations of for each model. 
In these figures we plot the boundary values of m° h for the parameters obtained as the 
solutions of our radiative symmetry breaking study. In each figure (a) the upper and 
lower boundaries of m° are drawn by using the all solutions obtained under the intial 
values shown in (24). In oder to show the ht(Mx) dependence of m' we classify the 
solutions into three classes and draw them separately. In figures (b) we plot the upper 
and lower boundaries of m" for the remaining solutions after imposing the additional 
condition m = 1 TeV. We also add the scatter plots of the all solutions corresponding 
to 2.4 TeV < u < 2.6 TeV in the same figures. They are represented by three kinds of 
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Fig. 4 Boundary values of m h a of the lightest neutral Higgs mass as a function of tan (3 in the model. 
Full data are used to draw 4(a). In 4(b) we impose rh(Mx) — 1 TeV. All solutions satisfying 2.4 TeV 
< u < 2.6 TeV are also plotted in 4(b). 

triangles corresponding to each h t (Mx). As a common feature in all models, we find that 
the larger h t (Mx) realizes the smaller tan (3 and then brings the larger contribution of 
the second term of Eq. (7). Thus the largest A(m t ) in the small tan/3 in Fig. 2(b) gives 
the largest m° h . Although X(m t ) in the extra U(l) models can be smaller than the one 
of the NMSSM as shown in Fig. 2(b), the boundary values of m° is larger in the extra 
U(l) models than in the NMSSM by a few to ten GeV. This is mainly due to the extra 
contribution to Eq. (17) coming from the singlet colored fields (gi, g~i). Since the existence 
of this contribution is the basic feature of the present extra U(l) models, the boundary 
value of is generally expected to be larger than the one of the NMSSM inspite of the 
running feature of the Yukawa coupling A. This one-loop effect is large enough to cancel 
the difference of X(m t ) in the second term of Eq. (7). In our studying parameters space 
the largest value of m° h is 



m 



o < 



156 GeV(NMSSM), m£ < 164 GeV(NMSSM), m£ < 158 GeV(NMSSM). (27) 



By comparing (a) and (b) in Figs. 3~5 we can get the tendency how the solutions are 
restricted when we reduce the parameter space. The change of m and u mainly affect the 
one loop contribution through the mass matrices (14) and (15). 
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Fig. 5 Boundary values of m h o of the lightest neutral Higgs mass as a function of tan (3 in the r\ model. 
Full data are used to draw 5(a). In 5(b) we impose rh(Mx) — 1 TeV. All solutions satisfying 2.4 TeV 
< u < 2.6 TeV are also plotted in 5(b). 

In Fig. 6 we plot the boundary value of m? h for u in the £_ model. This shows the 
tendency that the larger u gives the larger value of m|J. This is expected from the one- loop 
contribution of the extra singlet colored fields {gi,g~i)- From this figure we can read off 
the relation between mz> and m\ by using m%, ~ g\Q 2 s u 2 /2. The lower bound of mz> in 
Fig. 6 is about 600 GeV where we used gE{mt) — 0.36. The conditions (i) and (ii) also 
determine the lower bound of u in the extra U(l) models. 

Finally we give a few comments on some points related to the extra matters. We also 
studied the case (B) of the extra matter contents for the same parameter settings as the 
above study. In that case, as a common feature we can find, it becomes rather difficult 
to satisfy both of the radiative symmetry breaking conditions and the phenomenological 
conditions (i) to (iii) compared with the case (A). The number of solutions in the case 
(B) is drastically less than in the case (B). Since the value of g 3 (Mx) increases, h t {m t ) 
and k(m t ) becomes larger. In fact, the initial value of h t in the wide region such as 0.2 < 
ht(mt) < 0.9 results in only the small tan/3 (larger ht(mt)) solution such as tan/5 < 1.8. 
This also makes X(m t ) smaller. The larger tan/? solutions disappear and the value of \u\ 
is shifted upward. However, the upper boundary value of m? h behaves in the different way 
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Fig. 6 Boundary values of m h o of the lightest neutral Higgs mass as a function of u in the £_ model. 

between the NMSSM and the extra U(l) models. Although in both models m^o becomes 
smaller in the region of tan f3 >^ 2, the behavior is different at tan (3 <^ 2. In the NMSSM 
it is a little bit larger than the one of case (A). On the other hand, it becomes smaller than 
the one of case (A) by a several GeV in the extra U(l) models. Here we should remind 
the fact that even if \{m t ) is smaller m^o can be larger in the case that corresponding 
tan(3 is smaller. The difference in the RGE of A in both model is also important in this 
behavior. To have more confident quantitative results in this case we need to search the 
parameter space in the finer way. We also changed the number of gt) which couples 
to S in the superpotential (1) in the case (A). If we decrease this number from three to 
one, the boundary values of the allowed become larger. This reason is considered as 
follows. Although this decrease reduces the number of fields contributing to the one-loop 
effective potential, this also decreases the N g value in Eq. (10). As a result the larger k 
and A are realized at the low energy region. The larger k also brings the larger u. The 
contribution to the one-loop effect per a field can be larger. Thus the decrease of the 
number of (gi,g~i) which couples to S causes the increase of at not only the tree level 
but also the one-loop level. 

5. Summary 
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There are two well-known low energy candidates to solve the \i problem in the MSSM. 
These are the NMSSM and the extra U(l) models. We have estimated the upper bound 
of the lightest neutral Higgs mass in both models. Apart from a Higgs coupling \SHiH 2 , 
there is a typical coupling kS 3 in the NMSSM and kSgg in the extra U(l) models. In 
the NMSSM k plays a crucial role in the evolution of A which dominantly determines the 
tree level mass bound of the lightest neutral Higgs scalar and in the radiative symmetry 
breaking. In the extra U(l) models the introduction of the extra colored fields g,g and 
its coupling with the singlet Higgs S are crucial to cause the radiative symmetry breaking 
at the weak scale successfully. This coupling can also affect the running of the coupling 
constant A. We focussed our attention on these points and estimated the the upper bound 
of the lightest neutral Higgs mass in both models. In this estimation we additionally im- 
posed some phenomenological constraints related to A and the VEV of S coming from, 
for example, the mass bounds of the charginos, the charged Higgs scalars and the Z' 
boson. We solved the minimum conditions of the one-loop effective potential improved 
by the RGEs for the couplings and soft supersymmetry breaking parameters whose initial 
conditions are taken in the suitable region. We estimated the upper bound of the light- 
est neutral Higgs scalar for the parameters which bring the phenomenologically correct 
potential minimum. Its tree level contribution due to A can be smaller in the extra U(l) 
models than in the NMSSM. However, there is the extra one-loop contribution originated 
from the Yukawa coupling kSgg and this makes its upper bound larger in the extra U(l) 
than in the NMSSM by a few to ten GeV. It is interesting enough that the upper bound 
of the lightest neutral Higgs scalar in the extra U(l) models is not so different from the 
one of the NMSSM. The extra U(l) models may be an equal candidate to the NMSSM 
for the experimental Higgs search. 

This work has been partly supported by the a Grant-in-Aid for Scientific Research 
from the Ministry of Education, Science and Culture(#11640267 and #11127206). 
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